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ABSTRACT 

Among a number of gamma ray bursts whose host galaxies are known, GRB971214 stands out for its high 
redshift z > 3 and the Lya emission line having a P-Cygni type profile, which is interpreted to be a direct con- 
sequence of the expanding supershell. From a profile fitting analysis we estimate the expansion velocity of the 
supershell v exp = 1500 km s" 1 and the neutral column density Nhi = 10 20 cm" 2 . The redshift z = 3.418 of the host 
galaxy proposed by Kulkarni et al. (1998) has been revised to be z = 3.425 from our profile analysis. 

The observed Lya profile is fitted well by a Gaussian curve, which yields the Lya luminosity Li ya = (1.8 ± 
0.8) x 10 42 ergs s" 1 . Assuming that the photon source is a giant H II region, we deduce the electron number 
density in the H II region n e = (40 ± 10)(j^-) 5 ( 10 q pc )~ LS cm" 3 , which corresponds to the illumination by about 

10 4 05 stars. We estimate the star-formation rate to be Rsf = (7 ± 3)M Q yr" 1 with the internal and the Galactic 
extinction corrected. 

The theory on the evolution of supernova remnants is used to propose that the supershell is at the adiabatic 



E 53 y rs ' 



and the density of the ambient medium 



1/2 

phase, with its radius R = 18 E 5 j pc, its age t = 4.7 x 10 

ti\ =5.4 £53 1/2 cm" 3 , where £53 = Zs/10 53 ergs. And we estimate the kinetic energy of the supershell to be = 
7.3 x 10 52 £53 ergs. These values are consistent with the hypothesis that the supershell is the remnant of a gamma 
ray burst. 

We note similarities between supershells found in nearby galaxies and remote primeval galaxies, and propose 
that the gamma ray burst may have occured in a giant H II region whose environment is similar to in star forming 
galaxies. 

Subject headings: gamma rays: bursts — galaxies: individual (GRB971214) — galaxies: starburst — line: 
profiles 



1. INTRODUCTION 
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' Gamma ray bursts (hereafter GRBs), located at cosmologi- 
£>cal distances, form a group of the most luminous objects in the 
? '"jUniverse. A number of GRBs were observed with their host 
rNgalaxies, of which intensive spectroscopic and imaging obser- 
vations have been performed using the Hubble Space Telescope 
and the Keck telescopes. GRB971214 is one of such objects 
with a very high redshift z > 3 and also is worth a particular 
attention in following points. 

Firstly, even after the optical transient region had faded away, 
we can marginally detect a bright spot in the HST image, and 
the continuum and emission line fluxes from this spot over- 
whelm those from the remaining part of the host galaxy. Sec- 
ondly, the ultraviolet spectrum of GRB971214 illustrated in 
Fig. 1 shows a flat UV continuum which is often found in star 
forming galaxies. In addition, the Lya emission line has a black 
absorption trough in the red part of the emission peak. These 
facts imply that the UV spectrum is formed in a star forming 
region which is surrounded by a thick and expanding medium 
of neutral hydrogens. We note that the location of the GRB 
afterglow coincides with the star forming region or the bright 
spot, which leads to the proposal that the GRB occurs in a star 
forming region, and in this Letter we deduce the physical envi- 
ronment of the GRB from its spectrum. 

The P-Cygni type Lya emission in the spectrum of primeval 
galaxies has been often attributed to an absorption effect by a 
galaxy not associated with the primeval galaxy but interven- 
ing accidentally in the line of sight. It has been regarded as a 



damped Lya absorption that occurs in the vicinity of the source 
galaxy. In order to check this possibility, we calculate the prob- 
ability for observing an intervening galaxy in front of the GRB 
host galaxy, which is none other than the optical depth for see- 
ing a galaxy between the GRB host galaxy at z = 3.425 and the 
place that corresponds to v exp in the Hubble's expansion law. 
The optical depth is simply expressed by 



r = n„(\ +z) aL, 



(1) 



where the comoving volume number density of normal galaxies 
at z = 0, n g ~ 0.02/r 3 Mpc" 3 (Im 1995), and the path length L is 
estimated to be L = v exp /H with H being the Hubble constant at 
the redshift and given by 

:// [ftM(l+z) 3 + ftA] 1/2 , 



(2) 



where the cosmological density parameters are VLm = 
8ttGpo/3Hq, 51a = A/3// ( 2 , and the Hubble parameter //(> = 
65 km s" 1 
7r(riolO/i~ 
of 10 kpc. 
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Mpc" . Here, the cross section is given by a = 
1 kpc) 2 , where rio is the typical galaxy size in units 
A direct substitution yields the optical depth r 
0.0023r 2 for Q M = 1/3 and CL A = 2/3, and r = 0.0025r 
VLm = 0.3 and = 0. This indicates that if the average size 
of galaxies at z ~ 3 is not large, then it is highly improbable 
that the damped absorption in the spectrum of host galaxy of 
GRB971214 is formed by a galaxy intervening accidentally. 

This leaves us to consider an alternative hypothesis, accord- 
ing to which the P-Cygni type profile of Lya is formed by 
the the expanding supershell that surrounds the star forming 
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Environment of GRB97 1214 




FIG. 1 . — Profile of the Lya emission line and its P-Cygni type absorption (solid line) with the best fit profile (dotted line). The emission has a Gaussian profile 
whose width is 5.0 , line center flux /(A = 5380.8 ) = 0.675 fily. The line center frequency A = 5380.8 gives z = 3.425. The P-Cygni type absorption is best fitted by 
the column density of the supershell Nhi = 10 20 cirT 2 and the expansion velocity v exp = 1500 km s _1 . The dotted vertical line denoted by 'A' is the solid line center 
of our result, the line denoted by B is the location of the P-Cygni absorption trough, and the 'C vertical line is provided by Kulkarni et al. (1998). Comparing with 
Lya we also show SillAl 190 line, but the S/N ratio is not good. 



region in GRB971214 and is the remnant of the GRB prece- 
dent to GRB971214. In order to check this possibility, we 
now calculate a number of GRB event in a galaxy at z — 3.4. 
Assuming the supernova rate is proportional to the starform- 
ing rate, Sadat et al. (1998) calculated the supernova rate 
at z « 3.4, T SN ^ 0.011 x 10 6 h\ 5 SNe Myr" 1 Mpc" 3 , where 
Hq = 65 /i65 km s" 1 Mpc" 1 . Accepting the concept that the GRB 
rate traces the massive star formation rate, Woods & Loeb 
(1998) showed T CRB ~ 10 -6 r SA ,, where fl M = 0.3, n A = 0.7, 
and Ho = 65 km s -1 Mpc -1 . Therefore, T CRB (z ~ 3.4) ~ 
0.01 1 h\ 5 Myr -1 Mpc -3 . Adopting the number density of galax- 
ies at 2.0 < z < 3.5, $* = 1.76 x 10 -3 h\ 5 Mpc -3 (Pozzetti et 
al. 1998), we can get the GRB rate per a galaxy at z ~ 3.4, 
Tijrb — 6 Myr -1 . Therefore, the number of GRB events per a 
galaxy during 10 4-5 yrs is Ngrb = 0.06 <~ 0.6. Moreover, the 
beaming factor, if exists, can increase the event rate by another 
factor of ten and the number of GRB events per a galaxy dur- 
ing 10 4-5 yrs is Ngrb = 0.6 ~ 6, which makes our supershell 
hypothesis more probable and alternative suggestion. 

It is noticeable that the similar P-Cygni features are observed 
in primeval galaxies and nearby star forming galaxies. Lee 
and Ahn (1998) proposed that the features might be caused 
not by an overlapping intergalactic medium but by the expand- 
ing medium enveloping the star forming region in the galaxies. 
And this concept may be applied to such a remote star-forming 
galaxy as the host galaxy of GRB971214. However, we can not 
exclude the possibility that the multiple supernovae explosion 
may also result in the expanding shell. 

In this Letter, we adopt the hypothesis of GRB-driven super- 
shell and perform a profile fitting analysis to derive physical 
parameters characterizing the expanding supershell of neutral 
hydrogens surrounding the star forming region of the GRB host 
galaxy. 

2. IMAGES AND SPECTRUM OF GRB971214 

GRB971214 was detected at 9 UT, December 14, 1997 
(Heise et al. 1997), and its optical counterpart twelve hours 
after the burst (Halpern et al. 1998). With a total fluence 
of 1.09 x 10 -5 ergs cm -2 (Kippen et al. 1997) and the mea- 
sured redshift of z = 3.418 (Kulkarni et al. 1998), its energy 



release is estimated to be ~ 3 x 10 53 ergs in 7-rays alone un- 
der the assumption of isotropic emission, £Iq = 0.3, and Hq = 
65 km s -1 Mpc -1 . 

In Fig. 1 is shown its ultraviolet spectrum redshifted to the 
optical band and obtained by Kulkarni et al. (1998) using the 
Keck telescope. It is characterized by a flat UV continuum that 
is typically found in the spectra of star forming regions. It is 
also seen that the Lya emission has a P-Cygni type profile, 
which is frequently observed in the astronomical objects near 
or far (Lee & Ahn 1998). Lee and Ahn (1998) proposed that the 
P-Cygni type Lya line is formed when the Lya photons emitted 
in the central super star-cluster are radiatively transferred in a 
HI supershell that are optically thick and expanding. 

In this work, the photon source is assumed to be the H II 
region that may contain 10 4 O stars. According to Marlowe 
et al.(1995), nearby starbursting dwarfs are inferred to contain 
a similar number of OB stars, when considering Ha luminos- 
ity. So this can be thought to be neither entirely new nor ex- 
treme assumption for galaxies of higher redshifts. Furthermore, 
it appears less plausible that the Lya emission arises from the 
medium ionized by shocks produced by supernovae or hyper- 
novae. This is because the number density of the inner region is 
not sufficiently high to give the recombination time scale < 10 5 
years. 

3. IS THE PHOTON SOURCE SURROUNDED BY THE GRB 
REMNANT? 

3.1. Interpretation of the P-Cygni Absorption 

In this work we will consider the shell hypothesis, and derive 
the physical properties of the expanding neutral medium from 
the observed Lya absorption. 

We assume a Gaussian profile for the unobscured Lya emis- 
sion and convolve it with a Voigt function with the center dis- 
placed by the expanding velocity that will be determined by the 
fitting procedure. In principle, the effect of the frequency redis- 
tribution by back-scatterings should be considered. However, 
we neglect this effect in this paper, because the S/N ratio and 
the resolution of the spectrum are not sufficiently good. For the 
continuum level, the blue part of Lya, which is more prone to 
extinction, is extrapolated from the red portion of the spectrum 
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given by Kulkarni et al. (1998). They quote F v = \lA{v/v R )° ! 
nJy with a = -0.7 ± 0.2, where F v is the spectral density at fre- 
quency v and vr = 4.7 x 10 14 Hz, the central frequency of the 
R band. 

We show the result in Fig. 1, where the dotted line represents 
the best fit profile, the solid line the observed profile, and the 
horizontal solid line the continuum level. The best fit expansion 
velocity of the supershell relative to the H II region is deter- 
mined to be v exp = 1500 km s" 1 , and the best fit line center op- 
tical depth to = 6 x 10 6 , which corresponds to Nm = 10 20 cm" 2 . 
The best fit Lya profile has the width of a = 5 and the line 
center flux f(\ = 5280.8 ) = 0.675 /iJy, which gives the un- 
obscured flux to be 9.1 x 10" 18 erg cm" 2 s" 1 and the systemic 
redshift z= 3.425. 

This is slightly larger than the redshift proposed by Kulkarni 
et al. (1998), who may have overestimated the absorption in the 
blue part of the Lya. However, the absorption trough is suffi- 
ciently remote from the line center in the velocity space only 
to erode the extreme blue part of the Lya emission. Hence, we 
prefer the redshift of z = 3.425 of GRB971214 to the redshift of 
z = 3.418, and subsequently other physical parameters need to 
be revised. 

Assuming a standard Friedman cosmology with Hq = 
65 km s" 1 Mpc" 1 and Slo = 0.3, the luminosity distance di = 
9.7 x 10 28 cm. Considering the Galactic extinction, the un- 
obscured Lya flux is corrected to be F Lya = (1.5 ± 0.7) x 
10" 17 erg cm" 2 s" 1 , where the observational error given by 
Kulkarni et al. (1998) is introduced. Therefore, for 
the assumed cosmology, the Lya line luminosity L Lya = 
(1.8 ± 0.8) x 10 42 ergs _1 . If there is no internal extinc- 
tion in the interior of the Lya source, this corresponds 
to n e = (l^iO^X^) 05 ^)" 1 - 5 cm" 3 or n e = (40 ± 

10)( r ^)°- 5 ( T g^)" 1 ' 5 cm" 3 , of which the ionization can be 

maintained by ~ 10 4 05 stars as the ionizing source. 

From the Lya luminosity, we can estimate the star-formation 
rate (Thompson, Djorgovski, & Trauger 1995) to be Rsf = 
(7 ± 3)M Q yr" 1 , with both the internal and the Galactic extinc- 
tion being corrected. This is consistent with the star form- 
ing rate given by Kulkarni et al. (1998) as a lower limit, 
Rsf = 5.2 M Q yr" 1 which was obtained from the rest-frame con- 



tinuum luminosity at 1 , 500 . 

Using the revised redshift of the GRB host galaxy, we refine 
other absorption lines in the observed spectrum. In Fig. 2, we 
show the spectrum of the GRB host in the UV regime. It is 
seen that the revised wavelengths are in good agreement with 
the absorption features. 

3.2. Physical Configuration of the Supershell 

We consider the dynamical evolutionary model of the super- 
nova remnant to derive the physical quantities of the shell. Ac- 
cording to Woltjer (1972), the supernova remnant has four evo- 
lutionary phases, that is, the free expansion phase, the Sedov- 
Taylor or adiabatic phase, the snowplow or radiative phase, 
and finally the merging or dissipation phase (see also Reynolds 
1988). 

According to Woltjer(1972), the radiative phase begins 
roughly when the expansion velocity of the shell becomes 



v = 300(-^l T ) 2 /' 7 ( 



1 cm 3 1 53 ergs 



1 / 17 kms" 1 , 



(3) 



where n\ is the number density of the ambient medium and E 
is the initial explosion energy. Since the expansion velocity 
v = v exp = 1500 km s" 1 of the supershell exceeds the velocity in 
the radiative phase by a large margin, we propose that the super- 
shell is in the adiabatic phase, which is described by the Sedov 
solution. 

According to the Sedov solution in a uniform medium of 
number density n\ in which we have the relation n\ = 3N/R, 



R = Q.92{ ) 1/4 f 1/2 cm, 

m H N 

v = 0.37(-^— ) 1/4 f" 1/2 cms" 1 , 
Wl = 3 .2(^) 1 / 4 r 1 / 2 cm" 3 , 



(4) 
(5) 
(6) 



where v is the expansion velocity of the supershell, R the size 
of the supershell, E the initial explosion energy, N the column 
density of the supershell, m H the hydrogen mass, and t the age 
of the shell. 
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FIG. 2. — UV spectrum of the GRB971214 which is smoothed by a Gaussian kernel with its width a = 5 . Solid horizontal line stands for the continuum level, 
the solid vertical bars represent the line centers whose velocity components correspond to the Lya absorption. The refined absorption centers which correspond to 
the systemic redshift of the GRB host are presented by the dotted vertical lines. The dashed vertical bars stand for the absorption lines suggested by Kulkarni et al. 
(1998) The transitions for the absorption are also denoted. 
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Environment of GRB97 1214 



Using the values v = v exp = 1500 km s 1 and N = 10 20 cm 2 , 
we get 

f = 4.7x 10 3 (— ) 1/2 yrs, (7) 
N20 

/?= 18(^1)1/2 pCj (8) 

JV20 

ni =5.4(- 2 ^) 1 / 2 cm" 3 , (9) 
£53 

where N 2 o = N/10 20 cm -2 and £53 = £/10 53 ergs. 

In Fig. 3 are shown the size and age of the supershell as 
well as the volume number density of the ambient medium 
with £53 being a free parameter. For a range of input energy, 
0.01 < £53 < 100, we get the possible range of the other pa- 
rameters 0.53 < m < 53, 2 pc < R < 180 pc, and 5 x 10 2 yrs < 
t < 5 x 10 4 yrs. 

From these values we estimate the total kinetic energy of 
the expanding supershell given by = 2TrR 2 NHiniHV 2 xp = 7.3 x 
10 52 £53 ergs. It is also noticeable that this kinetic energy can be 
comparable to those of the galactic supershell including those 
in Our Galaxy, NGC 4631, and M101 (Heiles 1979, Rand & 
van der Hulst 1993, Wang 1999). Furthermore, the P-Cygni 
Lya lines of the primeval galaxies also show the similar en- 
ergy scale. Thus, we propose that the supershell is the rem- 
nant of a hypernova or a GRB that had exploded earlier than 
GRB971214. 

4. SUMMARY AND IMPLICATIONS 

We have studied on the formation of P-Cygni type Lya in 
the spectrum of GRB971214, and found that there are at least 
three components in the system, i.e. the parsec scale remnant 
of GRB971214 itself, a giant H II region, and a supershell sur- 
rounding it. 

The giant H II region from which Lya emission originates 
is photoionized by a super stellar cluster whose total ionizing 
photons correspond to those emitted from about 10 4 05 stars. 
The existence of the P-Cygni type absorption plausibly implies 
that there exists a supershell surrounding the H II region. By a 



profile fitting procedure, we found out that the shell is expand- 
ing with a velocity of v exp = 1500 km s" 1 and its neutral column 
density Nhi = 10 20 cm -2 . We also revised the redshift of the 
Lya emission source to be z = 3.425, and its unobscured Lya 
luminosity to be L Lya = (1 .8 ± 0.8) x 10 42 erg s" 1 , which gives a 
more reasonable star formation rate to be Rsf = (7 ± 3) M yr" 1 . 

We also applied the theory on the hydrodynamical evo- 
lution of supernova remnants to the supershell surrounding 
GRB971214. Assuming a reasonable scale of the initial ex- 
plosion energy of the supershell, we propose that the supershell 

1/2 

is at the adiabatic phase, with its radius R = 18 £53 pc, its 
age t = 4.7 x 10 3 E^ 2 yrs, and the number density of the am- 
bient medium n\ = 5.4 E 5 ^ 2 cm" 3 , where £53 = £/10 53 ergs. 
And we estimate the kinetic energy of the supershell to be 
E k = 7.3 x 10 52 £ 53 ergs. 

It is noticeable that there are many astronomical objects 
showing the similar characteristics. Using the X-ray data, 
Wang (1999) has discovered the candidates of GRB remnants 
in M101, of which two exhibit similar physical characteristics 
to those of the GRB remnant in GRB971214. With the advent 
of 10m class telescopes, a large number of primeval galaxies 
are observed by applying several methods including the Ly- 
man break method (Steidel 1996). About 50 percent of Lya 
emission lines in their spectra show P-Cygni type profiles. It 
is suggested that these profiles are formed in an expanding me- 
dia surrounding the star forming region and a case study was 
performed for DLA 2233+131 in detail (Lee & Ahn 1998). 

Recently one of the most debated suggestions is the hyper- 
nova conjecture, according to which gamma-ray bursts occur 
in star forming regions. Our results strongly favor this model 
and more concrete evidence is expected as the sample of GRB 
spectra showing Lya emission becomes statistically significant. 

The author thanks George Djorgovski and Shri Kulkarni for 
kindly providing the optical spectrum of GRB971214. He also 
thanks Bon-Chul Koo, Kee-Tae Kim, Hee-Won Lee, Hwang- 
Kyung Sung, In-Su Yi, and Hyung-Mok Lee for their invalu- 
able discussions. The author thanks to the anonymous referee 
for his/her fruitful comments and suggestions. 




FIG. 3. — Radius(R) in pc, age(( 4 ) in 10 4 yrs, and the volume number density of the ambient medium^) with the variation of the input energy (£ 53 ) in 10 53 ergs. 
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